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ABSTRACT 


An error probability analysis 1s performed for a conventional noncoherent M-ary 
orthogonal frequency-shift keying (MFSK) receiver employing fast frequency-hopped 
(FFH) spread spectrum waveforms transmitted over a frequency-nonselective, slowly 
fading Ricean channel with partial-band noise interference . Each diversity reception 1s 
assumed to fade independently. The partial-band interference is modeled as a Gaussian 
process. The effects of wideband thermal noise are also included. The energy per hop is 
held constant; thus, as diversity increases, energy per symbol increases. Previous 
analyses considered only constant energy per symbol systems, however, practical military 
systems are likely to employ fixed hop rates. There is some performance enhancement to 
be obtained from implementing diversity 1n a conventional FFH/MFSK system with fixed 
hop rates, but partial band interference still results in significant degradation. 

Additionally, the performance of this FFH receiver is investigated over the same 
channel in the presence of partial-band tone jamming without diversity for the case of 
binary frequency-shift keying (BFSK) when both the signal and the jammer can fade 
independently. Performance when only a single jamming tone per hop slot is allowed is 
compared to that obtained when two jamming tones per hop slot are possible. When the 
jamming signal experiences Rayleigh fading there is very little degradation of the 


jammer's effectiveness as compared to when the jamming signal is not affected by fading. 
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I. INTRODUCTION 


A fast frequency-hopping (FFH) communication system is a subset of spread 
spectrum communications that utilizes a bandwidth greatly exceeding that required for 
the information signal alone. Frequency-hopping spread spectrum is fundamentally 
different from direct sequence (DS) spread spectrum in the technique of signal generation 
and recovery. This thesis focuses on frequency-hopping spread spectrum systems 
because of their practical military importance. As spread spectrum systems grow more 
popular and occupy wider communications bandwidth, the likelihood of both hostile and 
non-hostile sources of narrowband interference also grows. It is important for mission 
planners to have a reliable estimate of the degradation that their communication systems 
will suffer under partial-band noise interference as well as tonal interference. Fast 
frequency-hopping spread spectrum techniques have evolved to counter the threat of 
intentional jamming [Ref. 1]. This thesis presents an error probability analysis of the 
conventional fast frequency-hopped orthogonal M-ary frequency-shift keying 
(FFH/MFSK) receiver with noncoherent detection for communications over channels 
with Ricean fading of the signal and partial-band noise interference. The influence of 
partial-band tonal interference on the binary orthogonal frequency-shift keying 
(FFH/BFSK) receiver without diversity where both the signal and jammer may fade 


independently is also considered. 


A. FFH/MFSK 


Frequency-shift keying is popular as a signaling scheme because it allows for 


noncoherent reception of the signal. A typical MFSK signal set can be expressed as 


si {2 A-cos(2nfot + 01) (1) 
so(t) = ¥2 A,cos [2n(fo + Af)t + 92] (2) 

and so on to 
su() = {2 A-cos[2n(fo + (M— 1 Aft Ou] (3) 


where f, is the lowest signal tone, Afis the tone spacing, and 9,, 1s the phase associated 


with each tone. This set is then modulated by a carrier that varies pseudorandomly in 
frequency for transmission as a frequency-hopped signal. 

All communication signals suffer from interference and noise. One method to 
overcome this degradation is to transmit the signal more than once thus providing a form 
of diversity. Fast frequency-hopping employs this redundancy as well as deliberately 
spreading the bandwidth. Frequency-hopping systems in which several symbols are 
transmitted per hop are considered slow frequency-hopping, while those transmitting 
several hops per symbol are considered fast frequency-hopping. Our FFH/MFSK 
transmitter performs L hops per data symbol, which results in a diversity of L. At the 
receiver the dehopped signals are recovered noncoherently by two correlators in phase 
quadrature with the dehopped signal waveform. The correlator outputs are sampled every 


T,, seconds, where 7, is the hop period. The sampled output of each correlator pair 1s 


squared and then these outputs are summed L times to obtain the decision statistic for 
each branch of the M-ary detector. The largest signal detected is selected as the 
transmitted symbol. A typical receiver structure is shown in Figure 1. 

Proper reception and demodulation of the spread signal depend on the recovery of 
several pieces of timing information. Both the sender and receiver need the same 
pseudorandom sequence operating synchronously. Also required are the symbol period 
and the hop period. In practice these are estimated from the received waveforms. This 


thesis assumes that this information is recovered without error. 
B. MULTI-PATH EFFECTS 


The losses experienced by the signal during propagation is worthy of an entire study 
in itself. However, it is useful to make some general observations on the composition of 
the received signal. It is possible and even likely that the received signal arrives at the 
receiver after transiting a variety of different paths. Signals traveling a longer distance 
arrive delayed relative to the direct path signal. This leads to multi-path effects. The 
magnitude of the multi-path effects depends on the magnitude of the delayed signal 
strength versus the direct path signal strength. It is common to consider the sum of all 


delayed signals as the diffuse component of the received signal. [Ref. 2] 


De-Hopped MFSK Input + Noise + Interference 
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S2quad(t) 





Sample every Thop 


Figure 1. Typical MFSK Receiver Structure 
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The consequence of multi-path reception is to cause the signal to fade in a time 
varying fashion. We can broadly characterize the channel conditions by examining the 
magnitude of the direct signal power to the diffuse signal power. The Ricean channel is 
the general case. Channels that have nearly all the received signal energy in the direct 
component, 1.e., direct-to-diffuse ratio greater than 100, have essentially no fading. In the 
limit, an infinite direct-to-diffuse ratio implies a Gaussian channel. Channels that have 
nearly all the received signal energy in the diffuse component, i.e., direct-to-diffuse ratio 
less than one, have strong fading. A direct-to-diffuse ratio of zero implies a Rayleigh 
channel. For direct-to-diffuse ratios between these extremes the channel experiences 
Ricean fading. This thesis examines the performance in each of these broad categories. 

In addition to these broad categories, the time varying nature of the channel can be 
described as slow or fast. In this thesis the channel properties are assumed to be constant 
over the duration of a hop and, therefore, slowly varying. Further, the channel may 
introduce some signal distortion arising from the treatment of sinusoids comprising the 
signal set within a hop differently. This distortion is characteristic of frequency selective 
channels. However, it is reasonable to assume that the signal sinusoids experience the 
Same multi-path effects. This is the case in frequency nonselective channels. One 
measure of this phenomenon is the coherence bandwidth of the signal. The coherence 
bandwidth is the frequency range over which the signal frequencies pass through without 


distortion. This can be summarized mathematically as 


(A). = 7 (4) 


where (Af), 1s the coherence bandwidth and 7,, is the multi-path spread of the channel. 
A frequency nonselective cannel displays a coherence bandwidth that is larger than the 


signal bandwidth. The rate of fading is related to 
l 
(A). =F (5) 


The coherence time 1s (Af), and B, is the Doppler spread of the channel. Slowly 
fading channels display a large coherence time or, conversely, a small Doppler spread. 
These descriptors are discussed further in the system description.[Ref. 2] The use of 
diversity to mitigate the multi-path effects for conventional MFSK has been widely 
investigated [Refs. 2, 3, 4]. 

When analyzing performance it is important to distinguish between a constant energy 
per hop system and a constant energy per symbol system. As diversity increases, the total 
symbol energy in a constant energy per hop system increases, while the symbol rate 
decreases. However, in the constant energy per symbol system, increasing diversity, L, 
implies decreasing hop duration, 7,.. Hence, a constant data rate 1s maintained, but the 
energy per hop is reduced. Since many practical military communication systems employ 
a fixed hop rate and a variable data rate, the constant energy per hop assumption is more 
logical [Ref. 5]. 


C. CONSTANT ENERGY PER HOP SYSTEMS 


Previous work has examined the performance of the noncoherent MFSK receiver in a 


Ricean fading channel with partial-band noise interference and constant energy per 


symbol [Ref. 3]. Additionally, the performance of the conventional BFSK receiver under 
partial-band noise jammin without the multi-path effect of fading is analyzed in [Ref. 6]. 
However, these investigations do not consider constant energy per hop signaling. The 
performance of several different types of diversity combining receivers, including the 
conventional receiver, utilizing constant energy per hop is simulated for Rayleigh fading 
channels in [Ref. 7]. 

Motivating this thesis is the uncertain degree of improved performance offered by 
more elaborate receiver designs, such as the noise normalized receiver, over the 
conventional receiver for variable data rate systems. [Refs. 3, 4] The expense and 
complexity of a more elaborate receiver may not be justified in some circumstances when 
utilizing a constant energy per hop system. 

The constant energy per hop scheme is of practical value because it allows the 
potential for an adaptive signaling scheme in which the sender and receiver can optimize 
the data transmission rate. In poor environments, diversity can be increased at the 
expense of lowered data rate. In favorable environments, the level of diversity, L , can be 
lowered to accommodate a higher data rate. These adaptations will ideally not require 
any hardware modifications and will be transparent to a channel observer. With this 
adaptive scheme in mind , this thesis examines the improvement offered by varying the 


level of diversity possible when the jammer is sub-optimal. 


D. TONE INTERFERENCE 


Another type of narrowband interference is tone interference. This can consist of a 
single interfering or multiple interfering tones. In this thesis, a performance analysis for a 
FFH/BFSK receiver without diversity over fading channel conditions similar to those 
assumed for noise interference is considered. Since the interfering tone or tones are 
signal-like in nature, they too can suffer multi-path effects. This analysis considers the 
effects of fading on both the signal and jammer. In a FFH/BFSK system an intelligent 
jammer can potentially cause more degradation by splitting his available power over both 
the signal tones. The degree of communication impairment of the single tone interference 
per hop versus two interference tones per hop strategy is also considered. Clearly, the 
greatest performance degradation occurs when the interfering tones correspond exactly to 
the various frequency-hopped symbol tones. Tone jamming where the tones do not 
correspond exactly to the various frequency-hopped symbol! tones are not considered in 


this thesis. 


Il. PARTIAL-BAND NOISE INTERFERENCE 


A. SYSTEM DESCRIPTION 


The partial-band noise interference, either intentional or unintentional, considered in 
this thesis is modeled as additive Gaussian noise and, when present, is assumed to be in 
all branches of the MFSK demodulator for any reception of the dehopped signal. 
Thermal noise and other wideband interferences which are also assumed to corrupt the 
signal are modeled as additive white Gaussian noise. Only the signal is assumed to be 
affected by channel fading. The smallest spacing between frequency hop slots is assumed 
larger than the coherence bandwidth of the channel, hence, each dehopped signal fades 
independently [Refs. 2,8]. As discussed in the Introduction, the signal bandwidth is 
assumed to be much smaller than the channel coherence bandwidth, and the channel 
coherence time is assumed to be much larger than the hop duration or, equivalently, the 
hop rate is assumed to be large compared to the Doppler spread of the channel. The first 
assumption implies that the channel is modeled as frequency-nonselective, while the 
second implies that the channel is slowly fading. The signal channel is modeled as a 
Ricean fading channel, hence, signal amplitude is a Ricean random variable [Refs. 2.8]. 
For Ricean fading, the total signal power consists of a direct signal component and a 
diffuse signal component, and the strength of the fading channel is characterized by the 


ratio of the direct signal component power to the diffuse signal component power. 


The symbol rate is R,. For MFSK with M order modulation, the corresponding bit 
rate is -%,= log ,(M). For Z hops per symbol, the hop rate is R,= LR,. The spread 


spectrum bandwidth, W, is considered very large compared to the hop rate. 


B. PARTIAL-BAND NOISE JAMMING ANALYSIS 


1. Problem Development 


The partial-band noise interference when present is assumed to be in all branches 
of the MFSK demodulator and affects each chip of the dehopped signal with probability y 


where y is the fraction of the spread bandwidth being jammed. Hence, the fraction of the 


spread bandwidth not being jammed is 1—y. If the average power spectral density of 
the interference is N,/2 over the entire spread bandwidth, then the power spectral 
density of the partial-band interference when present is N,/2y. The power spectral density 
of the thermal noise and other wideband interferences which are modeled as additive 
white Gaussian noise is N,/2. Consequently, total noise power spectral density is N,/2 in 


the absence of partial-band interference; otherwise, 
Nr a (Ne) 
a gg a 6 
ee oy) © 


is the total noise power spectral density when narrowband interference is present. 


If the equivalent noise bandwidth of each detector branch in the MFSK 
demodulator is B Hz, then the noise power received in a given hop is 0; = 07 =NoB 
with probability 1-y when no jamming is present. When jamming is present, the total 


noise power in a given hop kis 0, =07+07= (™ +No| B with a probability of y. 


We assume that each receiver hop slot has the same noise equivalent bandwidth. The 


noise equivalent bandwidth of the receiver investigated in this thesis is B=R,. 


2. Probability of Bit Error 


When partial-band interference is present the probability of symbol error for a 


MEFSK receiver is 


ee, 
ig -5/ . }ra - "Ps (7) 


— 


where ?,(i) represents the conditional probability of a symbol error given that i of L hops 


of a symbol are jammed. Since each signal branch of the receiver is symmetric with the 


other branches, we can determine the probability of a symbol error , P, , by considering 


the signal to be present only in branch one of the MFSK demodulator. The outputs of the 


other branches are assumed identical and independent (11d). 


by 


For orthogonal MFSK the probability of a bit error is related to the symbol error 


oe 
P» == 2(M — Dn: (8) 


The energy per bit as a function of the symbol energy and the modulation order is 


E. 
log,( 





(9) 


Ey = 
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3. Probability of Symbol Error Under Partial-band Noise Interference 


Assuming the signal is present in branch one of the MFSK demodulator allows 
us to write the probability of a symbol error based on the conditional probability density 
functions that ij of L hops are jammed.The conditional probability density function for the 


output of the branch containing the signalis fy, (x1 1) where X, is the random variable 


that represents the output of the signal branch. The conditional probability density 


functions for the non-signal branches are fy,,(Xmli),m =2,3,4...,M where theX,,'s are 


the identically distributed random variables that represent the output of the branches that 


do not contain the signal. [Ref. 3] The conditional probability of symbol error is 


P (i) = 1- [Pfcolal Pf ealidem | de (10) 


for all m+ 1. 


Since the partial band interference may or may not be present 1n a hop, we must 
be able to differentiate between the two possibilities. Let the subscript n= 1, 2 denote that 
hop & of a symbol has interference and has no interference, respectively. The diversity 
summer acts to add together all Z independent hops in its branch. Of those hops 7 of them 
are jammed; hence, we can express the random variable at the output of the diversity 


summer (10) as 


L i L 
Xm = 2, Xnh, =n + D, Xp eee (11) 
k=] 1 


k=L-1 





a. Probability Density Function for the Decision Variable X,, 
The probability density function of the independent, identically distributed 
(iid) random variables X,,, , m= 2,3,...,M that represent the demodulator branch outputs 


not containing the signal for hop k of a symbol is 





xu (mk) = +. exp wy wut) 6 6 = 2.3,...,.M (12) 
20; 20; 


where u(e)is the unit step function. [Ref. 9] 


The Laplace transform of (12) is 





| 
eS) = (13 
‘ Deeds) 


Since each hop is independent, from (11), we can express the Laplace transform of the 


conditional density function, fy,,(xmli) , as 
i 
Fy, pe | Fn, (sli) | x | Fn, (sli) | (14) 


Direct inversion of (14) leads to an infinite series of confluent hypergeometric functions, 
but this proved difficult to program and slow to execute. To evaluate (34) we first invert 

the individual portions consisting of either all hops jammed or all hops not jammed. The 
conditional density function for the decision variable is then obtained by convolution. 


The Laplace transform pairs are 


i 


| Fx, (sli) | > | | aN exp nae u(x) (15) 


2 
2044 2041 








and 


13 


L-i esi 
| Fn (sli) | = te Goan &%P exp =] u(x) (16) 


k2 
2 d 2 h . ry . ry e 
where 6;,° and 6,.° represent the noise power in a hop experiencing interference and a hop 


experiencing only thermal noise, respectively. To recover fy,,(xm) , the two m-Erlang 


random variables above are convolved . Hence, the inverse Laplace transform of (14) is 


obtained from a convolution of the individual right hand sides (RHS) of (15) and (16) as 


: oe L-i-1 
| I | i “ae m=) nO lg 17 
Sxl Ir | e Ga) hd 265 py eta “oe DOr ng 


Applying the binomial theorem to the integrand of (17) and defining 


i L-i 
_(.1_)\(._) 
ted ea (L—i-1)!G@-1)! (18) 








we can express (17) as 


ae ie ee | male is 
fin Emly = J Xu : k Jskcot IKexp Soe at (19) 


GF Ne, 


Interchanging the order of integration and summation and integrating term by term, we 


get [Ref. 10, equation 3.351.1] 


Las) sail oe Bee (— Dik+i=- Dia -l 


L-1-| 
: — _1\k L-1-I-ke 70 
Frbimli)= 2 2 ( k I en PGB) (ke i-1—Dya—By 


where for brevity 


a= — (21) 


and 
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p= (22) 


20%) 

The special cases of all hops jammed (i=Z) and no hops jammed (i=0) are 
obtained directly from (15) and (16) respectively, since the analysis leading to (20) is not 
valid in theses cases. The remaining task is to evaluate the cumulative distribution 
function of /,,,(x,,|2) that is required for computing the probability of symbol error, (10). 
The computation required for the non-signal branch contribution to the probability of 


symbol error is i" Ix,(%m). There are two avenues to compute this, the first based on 


direct integration or, equivalently, the second based on an inverse Laplace transform of 





Fy, (sli)/s. 
Directly integrating (20), valid when 0 <i< TZ, we have 
: _ ee a 
i Pin(eml idem = <r Fy Py CY) k 
ies 1 =e ENNIS Ie 
x(k+i-1)! a 7" 7 exp(—BXm) ee 
(a-By” (Bo Pa os iar be 


: k+i-] ( a 1)! 1 2 > L—2-m poem 
% (21 )@ep) one iti 2 (L-2-I-iar | | ©? 


This equation is simplified by expanding the binomial coefficient and canceling like 


terms to obtain 
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L-i-| ; ink 
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q - L-i-1-k Bey) > pei k+i-l (ao Nee 
| PC had p p! Ti aie yt 2 (k+i-1-D'\@-0)” a 





r ~ exp (arin) 3, a” | (24) 
q=0 p: 


when 0 <i< JL. For the case of all hops jammed, ij = L, the cumulative density function is 





L 
m = m l xl! »,§ 


2 
2044 


See eas uf) (25) 
2o7, | 20 K!\ 26%, 
and for the case of no hops jammed, j = 0, it is 


L 
m _ m 1 xh} 58 
Jy Saatedde= J, ea (L-1)! oi | fe 





k 
L-1 
=1l-—exp;—— >». 1 a (26) 
aaP K!\ 263, 


A final useful form is obtained by making the linear transformation of the random 


variable X,,. Recognizing that a = me and B = 5: then 


0] 05 
ea 2 O53 N 
Sie? 8 ZF 
203 ; 


Then making the change of variables y = Bx, 1n (24) we get 
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(28) 


which is useful for numerical integration purposes. Applying the same substitution in (25) 


and (26), we get 


L-} 
Fy(ylL) = 1 —exp -y} % 0)" (29) 
and 
IL-1 k 
Fy(yl0) = 1 -exp rap 12) (30) 





for the case of all hops jammed and no hops jammed, respectively. 

Previous efforts to determine the cumulative distribution for the non-signal 
branches, the second technique, is based on a numerical inverse Laplace transform of the 
function A x Fy,,(sli) [Refs. 3, 4]. Despite the finite sums of exponential-like terms in 


(28), this second method proved more accurate for large values of L . 


b. Probability Density Function for the Decision Variable X, 
Continuing with the assumption that the signal is present in branch one of the 


MFSK demodulator accompanied by Gaussian noise , the random variable_X, is described 


i 


by the conditional density function fy, (x1!) . Given a signal amplitude of J2 a, and 


before diversity combining, the probability density function of X,, is 


I 
xy, 1, 14%) = ae) 1 (31) 


2 
k 20; ; 


OX 


-(x14-+ 2a?) & 2x1, j 
———_ Ee eer) 


where /Jo(e)is the modified Bessel function of the first kind and order zero. The channel 


is assumed to have Ricean fading; therefore, a, is a Ricean random variable with a 
probability density function given by 


-(a? +a?) 


er fol 22) wax) (32) 


Fa, fan) = = eXp 


where o” is the average power of the direct signal component and 207 is the average 


power of the diffuse signal component [Refs. 2, 11]. The total average signal power of a 


hop k of a symbol is a? +207 and is assumed to remain constant from hop to hop. 
The random variable X\, when hop k of a symbol has interference and no 


interference is denoted by X1,, and X\x,, respectively. The conditioning is removed by 


integrating the product of (31) and (32) with respect to a, from zero to infinity with the 


result 


ae fy a 
ee Af Sua | po| Sa aa 
2(o?, 4 20?) 2 \ 04, +207 Ce Om 


for xix, 20 where es implies a jammed hop and oye implies no jamming. Since each 


hop is independent, from (11) we see that the Laplace transform of the conditional 
density function for the signal branch, fy, (x11) , is 


Fy,(sli) = (Fx, (sl ° [Fx (sly (34) 


The Laplace transform of the density function describing a single hop, fy,,(x1x) , is 


_ Br) eat 
Fy,,,(s) = (+ By) +B, exp| Be, (35) 


where 


ae 
_ 2(20? + o? | ” 


When(35) to the c, power and the inverse Laplace transform is taken. the result 1s 
equivalent to taking the c, fold convolution, denoted by @c,,. Letting c, =i and 


c, = L -i, so that elements subscripted with a 1 correspond to jammed hops while 


elements subscripted with a 2 correspond to hops experiencing only thermal noise, then 


(Can 1 il 
Bx,X Lk», 


| ee 
(2cp,at2) ) 


kin (Xk, ) | ae 


xexp [—Bx,(X1k, + Dena) a ee (26, gage x ju, Jax.) (37) 


where /,,-; represents the modified Bessel function of integer order c,-1. [Refs. 3, 11] 


Unlike the case of the non-signal branch analysis, no analytic solution for the final 
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convolution for the decision variable X, has been found. Symbolically, the probability 


density function for X, at the output of the diversity summer is described by 


@(L-i) 


fr eil)=[fru, Oru, ld] @ [Feu Crus!) | (38) 


C. NUMERICAL EVALUATION OF PARTIAL-BAND NOISE JAMMING 


For levels of diversity up to L= 20, the conditional probability of symbol error is first 
computed based on the signal energy-to-thermal noise power spectral density ratio, 
signal-to-jammer power spectral density ratio, signal fading direct-to-diffuse ratio, and 
fraction of spread spectrum bandwidth jammed. These input parameters form the basis 
for the probability of symbol error which 1s converted to probability of bit error using 


(8). To summarize 


) ; 
a; direct signal power/ho 

C= ree (39) 
Oj noise power/hop 


Iq2 diffuse signal power/hop 
Ck = aa = = es (40) 


noise power/hop 


Eee average energy per bit (41) 
N, _ thermal noise power spectral density 
En _ average energy per hop (42) 
Ny, interference power spectral density 
Gk Ot, 
Ry = ea = ples 43 
ae 562 (43) 


and 0.0<y<1.0. 
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The computation of the probability of symbol error is achieved for each L by 
weighting and summing the probabilities of symbol error conditioned on 7 of L hops 
jammed. A numeric integration is performed over the product of the signal path 


contribution and the non-signal path contribution (10). 
1. Non-Signal Path Contribution 


The cumulative distribution function is computed using the analytic expression of 
finite summations given by (28) for 0 <i< ZL for small ZL. However, it is computationally 


faster and results in less round off error for large values of L if a numerical inverse 


Laplace transform is performed on x Fy,,(sli). The special cases of all hops jammed 


(i=L) or no hops jammed (i=0) can be computed directly from (25) or (26), respectively. 
The convergence of the inverse Laplace transform is accelerated by the Euler 
transformation [Ref. 12] usually taking about 60 terms to reach relative errors on the 
order of 10°. However, certain combinations of input parameters results in slower 
convergence. The inverse Laplace transform algorithm is limited to not more than 1000 


iterations. 
2. Signal Path Contribution 


For the special cases either of all hops jammed or no hops jammed, the required 


probability function f;,(x\ 17) is given by (37). In the case of all hops jammed n=], 


c, = L, and c, =0, and 1n the case of no hops jammed n=2, c, = 0, and c, =L. In both of 
these cases, (10) can be solved analytically, but the results are so complex that numerical 


evaluation is easier and more straightforward. When # 0 and i # L,(38) is evaluated by 
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the numerical inversion of (34) after which the probability of symbol error, (10), is 
evaluated numerically. 
The partial-band jamming fraction that yields the worst-case performance for the 


conventional receiver is obtained experimentally by computing the probability of bit error 


En ag Eb 
N. and N; 


fading severity, results are obtained for several values of the direct-to-diffuse ratio. For 


as a function of y for fixed values of . To cover the broad range of channel 


weakly fading channels, exhibiting a strong direct signal component, a direct-to-diffuse 
ratio of 100 is used. This is essentially a Gaussian channel. For strongly fading channels, 
considered as nearly Rayleigh channels, a direct-to-diffuse ratio of one is used. Typical of 


a moderate Ricean fading channel, a direct-to-diffuse ratio of ten is used . 
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D. PARTIAL-BAND NOISE JAMMING RESULTS 


To determine the worst case partial-band jamming performance of the FFH/MFSK 
receiver the probability of bit error is computed as a function of y. The process is repeated 
for levels of diversity up to twenty, which represents the upper end for practical systems 
in today's technology. Figure 2 through Figure 9 display the probability of bit error for 
various levels of signal-to-noise ratio, total jammer power, and modulation order. The 
figures are calculated with the assumption that the signal direct-to-diffuse ratio (R,) is 
constant during a hop duration. Figure 2 demonstrates the performance for moderate 
signal-to-noise ratio with a near Gaussian channel. The signal-to-thermal noise ratio, 
E,/N,. iS the signal energy contained in one hop for a conventional BFSK receiver. The 


table below summarizes the cases displayed in Figure 2 through Figure 9. 


TABLE 1. PARAMETERS FOR WORST CASE eS 
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Figure 2. Worst Case Partial-band Jamming 
E,/N,=13.35dB, E,/N=3dB, R,=100 


The surprising result is that despite increased levels of diversity, the performance 
does not greatly improve. There will be more to say about this later, but the same trend is 


observed in Figure 3 where there is Ricean fading of the signal . 
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Figure 3. Worst Case Partial-band Jamming 
E,/N,=13.35dB, E,/N=3dB, R,=10 


Another, surprising result 1s that increasing the modulation order to four or eight 
does not provide much improvement. Figure 2 and Figure 3 display pessimistic 
performance from the communicator's viewpoint. However, they assume a high level of 
jammer power. We next investigate the performance when the jammer power is just one 
tenth of the communicator's power when L=1 ( Figure 4 and Figure 5). Since we are 
assuming 1000 hop slots, this is equivalent to assuming a signal-to-jammer power ratio of 


-20dB for BFSK with no frequency-hopping 
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Figure 4. Worst Case Partial-band Jamming 
E,/N,=13.35dB, E,/N=10dB, R,=100 
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Figure 5. Worst Case Partial-band Jamming 
E,/N,=13.35dB, E,/N=10dB, R,=10 
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A comparison of Figure 4 with Figure 5 illustrate the influences of Ricean fading 
on the probability of bit error when the communicator enjoys a significant power 
advantage over the jammer. Next we consider a situation where the communicator has a 
larger initial signal-to-thermal noise ratio, i.e. 16dB, to see if at higher levels the jammer 
could mitigate the effects of diversity with partial-band noise jamming. Figure 6 and 
Figure 7 display these results. 


Pb 
0.250 















0.200 





0.150 





0.050 —_—~ : ! —— : : 
123 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 
Number of Hops 


Figure 6. Worst Case Partial-band Jamming 
E,/N,=16dB, E,/N =3dB, R,=10 
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Figure 7. Worst Case Partial-band Noise Jamming 
E,/N,=16dB, E,/N=10dB, R,=10 


To examine the extremes of channel fading on the received signal, the case of 
near Rayleigh fading is examined with a direct-to-diffuse ratio of one. Since fading is 
severe, we would expect poor performance even without jamming, but with some 
improvement added by the increased diversity. The signal energy per hop is examined at 


20 dB (Figure 8) and at 10dB (Figure 9). 
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Figure 8. Worst Case Partial-band Jamming 
E,/N,=20dB, E,/N =3dB, R,=1 
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Figure 9. Worst Case Partial-band Jamming 
E,/N,=10dB, E,/Nj=0dB, R,=1 
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In both Figure 8 and Figure 9 the improvement offered by diversity initially 
lowers the probability of bit error, but not to the extent that it can overcome the errors 
introduced by partial-band jamming. It should be noted also that the high 
signal-to-thermal noise ratio is received after passing through the near Rayleigh channel. 
Therefore, it has a diffuse signal component equal to the direct signal component. This 
represents a very poor channel for the communicator that 1s further degraded by 
partial-band noise jamming. 

The conclusion to be drawn from the preceding series of performance plots is that 
for a wide range of signal-to-thermal noise level, jammer power, and channel conditions, 
the increase in diversity in the FFH/MFSK receiver alone is not enough to overcome the 
effects of partial-band noise jamming acting in concert with multipath effects. The 
critical parameter then is the fraction of the spread spectrum bandwidth, y, over which the 
jammer power is spread. The Figures 10 through 17 demonstrate the fraction of jammed 


bandwidth that provides the worst case probability of bit error. 
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Figure 10. Fraction of Jammed Bandwidth 
E,/N,=13.35dB, E,/N =3dB, R,=100 
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Figure 11. Fraction of Bandwidth Jammed 
E,/N,=13.35dB, E,/N=3dB, R,=10 
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Figure 12. Fraction of Bandwidth Jammed 
E,/N,=13.35dB, E,/N =10dB, R,=100 
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Figure 13. Fraction of Bandwidth Jammed 
E,/N,=13.35dB, E,/N=10dB, R,=10 


o2 

















123 45 6 7 8 Q9 1011 12 13 14 15 16 17 18 19 20 
Number of Hops 


Figure 14. Fraction of Bandwidth Jammed 
E,/N,=16dB, E,/N =3dB, R,=10 
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Figure 15. Fraction of Bandwidth Jammed 
E,/N,=16dB, E,/N=3dB, R,=10 
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Figure 16. Fraction of Bandwidth Jammed 
E,/N,=20dB, E,/N =3dB, R,=1 
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Figure 17. Fraction of Bandwidth Jammed 
E,/N,=10dB, E,/N =0dB, R,=1 
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The worst case values y are determined by computing the probability of bit error 
as a function of y. Figure 10 through Figure 17 show that , in general y, depends on the 
hop energy-to-thermal noise ratio (E,/N, ), the hop energy-to-jammer power spectral 
density (E,/N,), the level of diversity (L), the modulation order (/), and the channel 
fading direct-to-diffuse ratio. Further, the shape of the worst case y is generally the 
same, showing an initial drop as L increases with a gradually lowering slope. It is 
important to note that the lowest levels of y are still above 10° which represent only one 
jammed hop slot when there are 1000 hops in the frequency-hopping system. 

That the increase in diversity is unable to mitigate the partial band jamming is 
also indicated by looking at the trend of performance for fixed values of y. Figure 10 
through Figure 17 show that performance improves linearly as diversity increases. This 
is not unexpected since the constant energy per hop system increases the total energy per 
bit as diversity increases, however, the change in slope of the performance curves is 
indicative of the effect contributing to the worst case analysis. Figure 18 through 


Figure 23 display the performance for fixed values of y. 
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Figure 18. Performance for fixed values of gamma 
E,/N,=13.35dB, E,/N =10dB, R,=10 


Figure 18 displays the performance in a moderate Ricean fading channel where 


communicator enjoys a fair signal-to-jammer power advantage. The performance 


improves linearly as diversity increases. However, the slope of the performance curve is 


significantly reduced for decreasing values of y. Further, we can see that for a given L 


there corresponds a worst case value for y which becomes less sensitive as diversity 


increases. Figure 19 displays the performance when the signal enjoys less power 


advantage relative to the jammer. The trend 1s that for lower signal-to-jammer power 


ratio performance improves more slowly. 
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Figure 19. Performance for fixed values of gamma 
E,/N,=13.35dB, E,/N=3dB, R,=10 


The worst case fraction of bandwidth jammed 1s influenced by channel fading 
(Figure 10 through Figure 17). Likewise, the slope of the performance curves is 
influenced by channel fading . As the fading becomes more severe, 1.e., near Rayleigh, 


the performance improves more slowly as demonstrated in Figure 2] through Figure 23. 
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Figure 20. Performance for fixed values of gamma 
E,/N,=13.35dB, E,/N=10dB, Ry=1 
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Figure 21. Performance for fixed values of gamma 
E,/N,=13.35dB, E,/N=3dB, R,=1 
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Figure 22. Performance for fixed values of gamma 
E,/N,=13.35dB, E,/N=10dB, R,=100 
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Figure 23. Performance for fixed values of gamma 


E,/N,=13.35dB, E,/N=3dB, R,=100 
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Ill. MULTI-TONE INTERFERENCE 


A. SYSTEM DESCRIPTION 


In addition to partial-band noise jamming, spread spectrum communications can 
experience tonal or narrowband interference in one or more hop slots. We can similarly 
expect that as spread spectrum bandwidth grows the likelihood of tonal interference will 
also grow. This interference may arise from hostile sources as well as non-hostile ground 
communications, satellite transponders, and radars. The multi-tone interference problem 
is similar to the partial-band noise problem. 

This chapter analyzes the probability of bit error arising from an "intelligent jammer" 
who knows the spread spectrum bandwidth, W, and the modulation order. This 
intelligent jammer is able to place a single interfering tone in one or more hop slots, but 
has no knowledge of the pseudonoise (PN) sequence driving the hopping pattern. This 
represents a very intelligent jammer. We then relax this situation and compare it with the 
performance of the jammer who places an interfering tone in none, one, or both of the 
symbol tone locations in a particular hop slot. 

It is likely that the interfering tones may experience multi-path effects in a manner 
similar to the signal tone. However, since the interfering tones and signal tone can arrive 
from different paths, the interfering tones will in general arrive at the receiver with a 
different direct-to-diffuse power ratio from that of the signal. In this thesis, the 


performance when both the signal and the jammer experience fading is evaluated. 
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B. MULTI-TONE INTERFERENCE ANALYSIS 


1. Problem Development 


First we consider the case where the jammer places at most one tone in each fast 
frequency-hop slot. Only the worst case situation where the jamming tones exactly 
coincide with frequency-hopped symbol tones is considered. The situation may be 


visualized with the aid Figure 24. 


FFH/BFSK 
N slots 





| tones jammed | 


Figure 24. Multi-tone interference in FFH/BFSK 
The FFH/BFSK receiver is a special case of the FFH/MFSK receiver described 


earlier. The orthogonal design 1s selected such that .#, 1s the system bit rate and each 
signal tone 1s placed -#, apart. In Figure 24 the bandwidth of a single hop, B, is 2.7, . 
The spread spectrum encompasses W Hz and is divided into N hop slots. Although 


numbered consecutively for illustration, the pseudonoise modulation will create a 
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random-like sequence to the hopping pattern. In this thesis, N=1000 is choosen as a 
representative number of hop slots. Further, the number of interfering tones, g, must be 


an integer 
sire Per (44) 


Since the jammer spreads his available power, J, equally among g tones, each 
interfering tone has a power of J, =J/q. Earlier work focuses on the simplifying case 
of no thermal noise [Ref. 13]. In the no thermal noise case an error can only occur if the 
hop is jammed and the the interfering tone power is larger than the signal power. In this 
instance the jammer can make most efficient use of his power by choosing the number of 


tones to jam as 


_ cle 
qg =INTG Z| (45) 


where P.. is the carrier power of the desired symbol and INTG[ ] represents the greatest 
integer less than the argument. [Ref. 13, pp.596-598] In this thesis, the probability of bit 
error is obtained using fixed values of g to determine a worst case scenario which 
includes the effects of both signal and jammer fading as well as thermal noise. The 
approach to this problem is to first detemine the probability of bit error for a conventional 


BFSK system without frequency-hopping and then generalize those results to the 


FFH/BFSK system. 
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2. BFSK Analysis with Single Tone Interference 


For the binary case we begin with the familiar expression for the probability of 


bit error for the noncoherent FSK receiver 


| etd eae car UNG > Rylsy]+Ps, x P,[Ry > Ros] (46) 


where R, and R, represent the demodulator outputs of the branch detecting signal s, and 


S, , respectively. Making the assumption of equally likely signaling, we have 


P,, =P, =+ (47) 
2 
and 
ae S1P,(Rs Sf OG Ss real (48) 
The received energy consists of signal energy, thermal noise, and, if jammed, the 
interference signal 
rO=s(H+nO+ssO, (49) 


where 5; 1s the ith signal tone, s, is the interference tone, and n(t) is additive white 
Gaussian noise with power spectral density N,/2. Since the receiver structure is 
symmetric, we can simplify the analysis problem by considering the case where branch 
two 1s jammed. 

First we examine the case where the signal and interference occupy different 


frequencies. Then 


P,(errorl|s,) = P,(R2 > R1\s1) (50) 
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We express the received waveform as 


r(f) = {2 a-cos (11+) + 2 aycos (wt + 0y) +n(t) 


(51) 


where a, and a, are modeled as Ricean random variables. As previously discussed the 


received amplitudes and phases of each tone is a random variable due to multipath 


effects. For the noncoherent detector design, the output of each branch can be expressed 


as anon-central Chi square random variable [Ref. 2]. For the signal branch 


(>, +2a? a 


Dot 


jo E8 Jry Jur 


l 
FrGilassii— 552 EXP = 


and for branch two containing the interference tone 


(+. +2037?) 


20? 


Pr (r21a9, 51) = exp _ 
20 


where o” = N.7,. The probability of error (50) is 


P,(R2 > R, lac, ayS1) = i IF PRR (rir21ay.a¢.8))dr2dr 


This represents the probability conditioned on the amplitudes of a, anda,. The 


orthogonal receiver structure and independent fading of the interference tone and the 


signal tone allow the joint probability density function to be separated. Hence. (54) 


simplifies to 


P(R2>Ri Sy [fac lacys1)| Je fe (r2lay.si dre fdr 
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(52) 


(53) 


(54) 


(55) 


The evaluation of (55) 1s simplified by applying the quadratic transformation to convert 


(52) and (53) to Ricean random variables [Ref. 11]. Ifr; =v7, then a =2v,. Using 
Vv) 


this substitution in (52), we obtain 


2 272 
Vv] (vi +2a27}) 2 acT pv) 
fiona ot) = ee p> /)| —-————- 


i * poo (56) 


Similarly for the branch output containing the interference tone 


2 2 2 
V2 (3 = 24573} WO fae 
fv,(v2lay,51) = 5 © AP — a Io a, u(v2) (57) 
20 20 0) 
Now (55) can be expressed as 
P,(V2 > Vils\) = [fv (Vv) lae,51)) I frs(v2las.si)av2 [av (58) 


The inner integral equation has no simple solution, but 1s represented symbolically as 


Marcum's Q function which 1s defined as [Ref. 14] 





O(a, B) = if vexp 7 + a)| Iofowr}as (59) 


Hence, (58) 1s expressed as 


Aaa: T.. 2 acTr 
P.(v2>vilsi)= [° O( 2 asTolo. 24) x *bexp| R2ETE|,| BacTon |, (60) 


2a- 


which can be evaluated to obtain [Ref. 15] 





l ole ayl inega! 
P(v2 > vy 151) = aL ag 2 = ; +o Sef, 27a) (oy 
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Equation (58) represents the contribution to the probability of bit error when the signal 
and jammer occupy different tones within a hop. It remains to determine the contribution 


due to a collocation of interfering and signal tones. That is 


P,(error |s2) = P,(R, > R2\|s2) (62) 


where 


Jins(r1) = 15 exp | 4 utr) (63) 


Since the signal and the interference are in the same tone location, the output of 


the in-phase and quadrature portions of branch two due to signal and jamming tones is 


{2 Ts[a-cos (8) + aycos (0)) (64) 
and 


2 Tslacsin(®) + aysin(6,)] (65) 


respectively. The phase of the jammer is generally different from that of the signal. 


Squaring and summing to form the branch output produces 


R20 — 04) = 27;| a2 +47 + 2aca,(cosOcos 0, + sinO sin 0, | (66) 


which is simplified to 


R3(d) = Didar + as + 2a,ayCOs 8 | (67) 


where 5= 9 -0,. Now, replacing 2a,;7,° with (67) in (53) we get 


1 (+. +273(a? +03) | jf = ae, 
fr(12\ac,.4),8,52) = 55 EXP a Io ae ana? 24.03) u(r2) (68) 


Now, 
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Feral ee R> Is2, 5) = |, fex(ralacsay, Ons) ) x iP PRM Is2)dr} lar: (69) 


Thus, the probability of error arising from collocation of the jammer and the signal is 
computed from a probability density function conditional on 6. In this case the 
non-signal branch is unaffected by the jamming tone. Substituting (63) and (68) into 


(69), we get 


7 
Jer O ie > R> Iso, 0) = ; xX exp F(a +a te 2a-asCOS 5) (70) 


This can be simplified into familiar terms by replacing 


14 


£ 
ao oS = 1 (71) 


and similarly for the jammer power term 
fe ma any — | (72) 
Then the conditional expression in (70) then becomes 
P(R\ > Ralsi,8)=(4Jexp{S(y, +yst+2/yeys cos6) | (73) 


The nuisance parameter 6 may be removed by multiplying by its probability density 


function and integrating the conditional density. That is 


PaCS ls) ) = _ ig exp | (y. +¥t+2 /yovs css] | dé (74) 
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where 6 is modeled as a uniform random variable. Since the integral is only over 6, then 


[Ref. 10] 





P(Ri > Ralsi) = exp {st +7] x be: i exp | Iyeyy cos} dd (75) 


which can be evaluated to obtain 





P,(Ry > RyIs1) = (1) x CXp -5(¥e +y)| x Io| [vers | (76) 


Therefore, the phase relationship between the signal and the jammer plays a significant 
role as evidenced by the modified Bessel function term. This makes physical sense. 
When the phase 1s the same or nearly so there is reinforcement of the signal tone by the 
interfering tone, but when the phase difference approaches 180 degrees there is 
destructive interference. 

This allows us to compute the probability of bit error in the absence of fading by 


substituting (61) and (76) into (48). 


Palaa= of f222, es) shal fe 
2 
sexp |e +a) |r) $22 (77) 


An equivalent expression using (71) and (72) ts 


Poy) = 41 -O( 7. J) +O( Jw. Jr) +exp 








> (y +vy}ol fies || (78) 





Using a QO function property [Ref. 15, pp. 396], we can reduce this to an expression 


containing just one Q function as 
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P gett) = 1) xo jasT jacTs | 
—| I acasTs —1]xex “Ta (42 + a3) (79) 
"fg? a 


Now, the amplitudes of the received interference tone and signal tone are Ricean random 


variables where 








2 
a;,+o 

fa,(as) = a | eng x by 2 la (80) 

Oy 207 O7 

Ac at +a2 Ace 

fa.(@c) = = exp }-| Sa | ¢ < Lo] > [ul@c) (81) 

O< ZG - Oc 

Where in each case the second moment can be expressed as 

a? = a2 +20? for i= Jie. (82) 


We first remove the conditioning on a,in (79). The technique is to take each piece of 
(79), multiply by (81), and integrate over all possible values of a,. Hence, for the first 


piece 





I of, au ? aa Vaasa a 
Ir of fat sc 7 = Ez <7} = j x | 22 laa (83) 
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2-2 2 m2 2 m2 

= Ooi; Te els ajTl, 

=a a. | ee 
627; +07 Gel ee Weis 


an | of | * 


bi 272 2 272 De 272 2 
O-ip,pto O-i, to O-ip,pts 


The integral of the product of the second piece and (81) is 


ee ( tenes aoe, a Cs aco 
exp ) -—~| az +. a4} po) ——— | x exp 4—-| —— | plo] —<lda. (85 
. n| 207 i) y 202 Ge ‘i Digs 62 ce 
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Making the substitutions 


=a aoe ; 
267 2G. Ga 





eal i aT; 
é =——) andl = j-, then 
Oc 





ye) 71 a ne 44 
B =e = ajT,o% +A60 . 26 Oc (86) 
4G 40%o% T,02 +07 
BY _ ja-ayT; (87) 
20 T;02 +67 
in (85), we can use [Ref. 10, equation 6.633.4] to obtain 
e -7(a2 +02) Pact, 
aces (lic. > eee ———————— (88) 
T,0.+60 2T,02 +26 T,02 +07 


for the second piece of (79). The remaining piece to be determined is 


lis 2 2 Ac a2 + a2 aAcAc 
[, exp |i (ai +2) } 4s exp) — 262 e 1, a da (89) 
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Using [Ref. 10, equation 6.631.4], we evaluate (89) to obtain 





4 272 272 
yearn 26 2( T20? + s?) 


Finally, combining (84), (88), and (90) with (79), we get the probability of bit error when 


the signal experiences Ricean fading as 


2 2 2 2 
ee Sols ro eT sy | aTs i 
2) oe hee . wee kor 


2( 02 eae 


oe O aiTs 27s i 
ea 8) | aie a ea pa 
ae ee oti Boe 





2(0273 +02] 
2 oi 
= exp 4—- 2 x 
4( 7302 + o?) 2( 720? + 3?) 
2 2 2 m2 2 
Ol Oouasg 
2| iB exp) Jt b x Io) > (91) 
2( T202 + 3?) Ge ee 


which 1s now conditional on a,. Expressing (86) in terms of 








Coe a2T; 
GE ONG a 
202T+ 2627, 


and the direct-to-diffuse power ratio 





(94) 


we get 





1 205 De, _ [2asTi/o? 
P,(a ee E, =n d Ep +2 [co |} 
o| [2asTalo? | 20» <P Cp g 
ae E,+2 Ep+2 Es (€5 +2) 
Ah Dee 
nf eof 8B] gy 


2ay Tylor 2p 2C ae lan in ee 
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| | -] “eu adie 2) 0 
Xx £0 
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Using [Ref. 15] 














in (95), we obtain 


4 265 | 2ajT;/0? i {|=} 
ray=Hi-gf 2D 2d Wael: |Iachges Oe 


2ajT 
en {-aT;,/20? } | +e] = =aT . oe al (97) 


oO *(E, +2) o(E, +2) 
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We now employ a similar approach to remove the conditioning ona,. Integrating 


the product of (80) and (97) over all possible receiver jammer amplitudes, a,, we get 


_l},_9 [7 J 
er pee es mn a 


saliasia Dye aimee bea }* : 2S, = | . 
where 
= seth (99) 
p= Pet _ Doth (100) 
a oF os ave 


An alternate expression for the probability of bit error in the case of both signal and 


jammer fading in terms of weighted QO functions is 


ogi l = Sy oe, 
reWi (fs /5,)| 
Ze; 206 op +Cy $3 aie, 
fe cl tof 25] 0 


However, since there is no closed form solution for the Q function, in practice it is easier 


to compute the performance using (98). 
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a. Special Cases of BFSK with Single Tone Interference 
The final aspect of the analysis is to consider how the derived result 
simplifies for special cases, both as a check on earlier work and to see limiting cases. As 
check on (98), when there is no signal the probability of bit error reduces to one half as 
expected. 
Next consider the case of no jamming. In this case, &;= C, = 0, and the 


probability of bit error reduces to 


—f 7 eS: 
r-(-45] exp} | (103) 


which is the usual BFSK result. 








Another limiting case is the performance when both the signal and the 
jammer experience Rayleigh fading. In this case both ¢, and ¢, are zero and all the 


received energy is in the respective diffuse component. In this case (98) reduces to 


_1f_6&+2 
p= 2 od 


This displays an inverse linear relationship similar to Rayleigh fading without jamming 
[Ref. 2]. Performance is improved if we change the case slightly to that where the 


jammer suffers Rayleigh fading while the signal experiences Ricean fading. Then 


] ( Eyt+2 os 
Pe oe E_pteyt2 exp] ( 
Here any direct signal component that reaches the receiver will serve to drive the 


probability of bit error lower in an exponential fashion. 
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Perhaps the most optimistic performance from the communicator's point of 
view is when the signal has no fading, but the jammer suffers Rayleigh fading. In this 


case, C,;=¢€,=0, and (98) reduces to 
5 =a 
=—1 106 
Ps ¢ xp | =| (106) 


On the other hand , the most pessimistic performance is likely to occur when 





the signal suffers Rayleigh fading, but the jammer has no fading. In this case, 


¢, =&,=0, and (98) reduces to 


| oo =6s 
rad tow { =} (107) 


This performance will be very poor unless there is a relatively large amount of diffuse 








signal energy received. 
3. FFH/BFSK Analysis with Multi-Tone Interference 
We now turn our attention to the application of (98) to spread spectrum and the 
partial band multi-tone jamming scenario. Suppose N is the number of hop slots in the 


FFH/BFSK system and the jammer transmits qg interfering tones. When we credit the 


jammer with only placing at most one tone per hop slot then 


v=(4) (108) 


is the fraction of spread spectrum bandwidth jammed where y is restricted to 


<y<1.0 (109) 


ale 


Now, 


ae 


. N- 
P,= (4) P,ferror|hop jammed] + =4 } p,ferrorlhop not jammed] 





where 


P,ferrorlhop not jammed] = P,[error|BFSK] = ts) exp i=} (110) 
b b 





To apply our earlier result for the BFSK receiver we must define our jammer parameters 
on a per hop basis. If 7, is the average total jammer power, then the average total 


jammer power per hop is 

Sb So , (111) 
which implies €, and C, are replaced by NE,/q and NGj/q , respectively, in (98). Now (110) 
is given as 


Py= (2 ee es 

2N E,+ RE +2 VEgt FE, +2 
7) (Brt2) [Gee te) | | 2 feeke, 
(Wf: ) ph clea an x Io NSEC Sakina 
2N 2(e,+%e) +2) Co aed tee Ch geome 








oe (Jen (5) a 


In the case of FFH/BFSK with one hop per bit, no fading, and no thermal noise 


(N,=0), an analytic expression for the worst case number of jamming tones, q,, 1S 
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(NYa\ _ 
1G | ag eek ad (113) 


One aim of this thesis is to determine if this limited analytic result for g, remains true for 
FFH/BFSK over Ricean fading channels in the presence of thermal noise. This is 


accomplished by comparing the performance based on q, and some offset values of g. 
4. FFH/BFSK Analysis Allowing Two Interference Tones 


We now extend the previous results to a situation which includes the possibility 
that both receiver branches contain an interfering tone. This 1s a relaxation of the 
intelligent jammer scenario employed thus far. In this situation, the probability of bit 


error may be expressed as 


P» = P,[{no hops jammed] x P,[errorlno hops jammed] 
+P,[one tone jammed] x P,[errorlone tone jammed] 


+P,[two tones jammed] x P,[errorltwo tones jammed] (114) 


The first two conditional probabilities were previously determined in Section 2. 
However, the probability of those events occurring is a sample without replacement 


situation. That is 


P,[one tone jammed] = P,[one tone jammed one tone not jammed] 


P,[0 tone jammed I1 tone not jammed] x P,[1 tone jammed |0 tone not jammed] 
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~ \QN7 0 N= 1) 7 \ON/\ON-1) ie) 
Similarly 
, q q 
ripe met} = (1-6) (1-5) 
[no tones jammed] aN aN oI 
PONE (ENO 
MN) ay Se Oe 
And lastly, the probability that both tones are jammed is 
P,[ both jammed] = P,[1 jammed 1 0 jammed] 
mec yf Ga 
ar \2n-1) mm 


All that remains 1s to determine the probability of error given that both branches are 
jammed. 
a. Extension of Single Interfering Tone Results 
Now that both branches of the BFSK receiver contain jamming tones, 
symmetry is restored . [Ref. 11] Hence, we can determine the probability of bit error 
supposing that signal s, is sent. An error occurs when the output of branch two 1s greater 


than the output of branch one and (54) yields 


P,[errorltwo tones jammed] = (1) x 


ab - (118) 


Tp Jaz + a + 2a-ascos ae: fis 10% (azt a; +2a-a,cos& 
(= ¢, | oe ee a 
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Alternatively, this may be expressed as 


2 
aT, 154 4¢ + a) + 2a-ajcosd 


Oe: Oo 





P,ferrorltwo tones jammed] = Q 


(1) 73( a2 +2a3 +2a,a,cos6) Tea; la? +a? +2a,a,cos 6 
7 exp )— 


x Io 
207 o* 


(119) 


The conditioning on the received amplitudes of the signal and the jammer is removed by 
multiplying (119) with (80) and (81) and integrating over all values of a, and a,. In 
addition, the conditioning on the signal-to-jammer phase difference, 6, must also be 


removed. 
In order to make (119) valid for FFH, we substitute (111) into (119). Now, 


symbolically the probability of error in the case of both tones jammed, P,, , 1s 


Pa=(1L) PP rtar taao 
2 20 pee a2+ Maj +2aca, ot cos 8 2 


lie (a2 + at+ 2a-a; cos 3) 
(!) 4) a 2 
20 


Teas zi 
; Thay S a2 + Na} +2aca, ea cosé | ¢da.dasdé (120) 
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Substituting our previous results for BFSK and (120) into (114), we get the 


total probability of bit error 


r= (Sa H(- LS) exp | | 


Ly ew 260 456s 
(2) =F 4 Ey t+ ey +2 = 








(36,42) (60+ 7) 2 o0FGs 
——— ERP ye ee 
tree ae Eg + SEs+2 Cp + Gye 


(4) or “| Pra(Go.Bs.bu CJ; ()) GAs 


Further simplification results by using the substitutions 





Ba (122) 
which implies 
dy, = “Se bda, (123) 
and 
2 
ail 
ie (124) 
which implies 
dy) = “bday, (125) 


Then using (122) through (125) in (80) and (81) 
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Ye* Sel, YeGe 
sSlai[sfE] om 





Sre(Ve) = A) exp 


and 


foo-()ool-22) iff] aan 


b. Special Cases when Allowing Two Interference Tones 
As a consequence of the triple nested numerical integration required by (120) 
in order to evaluate (121), obtaining numerical results when two tones can jam a single 
hop slot is computationally intensive. Based on the results when thermal noise 1s 
neglected, we expect the performance will be better than when only a single interfering 
tone 1s allowed. 
When there is no signal fading 
Th =p = Gp sel on =U) (128) 
and when there is no jammer fading 
Vi=y7— cy and e7— 0 G29) 
In the case of Rayleigh fading of the signal 
Ys =§» and C, = 0 (130) 


and 


when the jammer is Rayleigh faded. Two limiting cases allow evaluation of the most 


pessimistic and most optimistic performance. Optimistically, from the communicator's 
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point of view, when there is no signal fading and Rayleigh fading of the jammer, the 


probability of bit error reduces to 


(2N-q\(2N-q-1\__ {Ce (2N-q ) —Op 
P= (aN J IN—4 Jexp|*] +(0 ieee rex > 








2 


2N- 
(at me le q \ 2}(L) in [ fica ~ {o[ JB Yo + 2Nys/q +2 2Nysyulq cos8 | 
(2) exp So, +4Nyj/q+2,/2Nysys/q cos 5) 
x al (2Nyu/q Jo + 2Nys/q+2J2Nyevs/q cosd |e (132) 


which requires two numerical integrations rather than three. 
The most pessimistic performance result from Rayleigh fading of the signal 


and no fading of the jammer. In this case, 





ONE G (BN =e) Ce ONE =e 
Carre mor ols ]+(ah pena a a a ae | 
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Further simplification when neither the signal nor the jammer fade. In this 


instance (132) reduces to 


Poe Se) anet JE* 


(£)(PE8)[1-o( (2) 
cn | Lea +2nss nf JFREEM 
CIE) & {0 ner. fre +20rta 2 fire oo 
(Sex |(f) 4m +2 aR x) 


where a single finite numerical integration over the phase angle 1s required. In this thesis 
only (134) is evaluated. Further work in this area could consider the effects of Ricean 


fading in the both tones jammed case based on (120). 
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C. MULTI-TONE INTERFERENCE NUMERICAL PROCEDURE 


The crux of the numerical computation for this problem is the accurate and efficient 
computation of Marcum's Q function. Originally the O function described a shorthand 
notation for the probability integral representing the output of a correlation detector 
containing a radar target in the presence of narrowband Gaussian noise [Ref. 14]. 

The Q function, described by (59), has no closed form solution expressible in a finite 
number of terms. Frequently cited equivalent expressions containing an infinite series of 
Bessel functions, while valuable as analysis tools, do not ease the computational burden. 
The strategy is to test the input arguments for their magnitude, difference, and the 
presence of zeros. Based on the input arguments the program computes the value of the 
O function in an appropriate subroutine and returns to the calling program. 

First, the presence of a zero passed as an input argument is tested. If present the 
following simplification results 

O(a, 0) = 1.0 (135) 


and 


(0, B) = exp {-(£) 8} (136) 


The second consideration is the magnitude and difference of the input arguments. If the 
product of the arguments is greater than ten and their difference is greater than five, then 
the QO function is computed using an asymptotic polynomial. If the magnitude of the 
product is greater than 1000, then the Q function is computed by directly integrating (59) 


using a large argument approximation for the Bessel function. For the semi-infinite 
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integration required here, the Romberg technique is very efficient. The Romberg 
technique is a recursive implementation of Richardson's extrapolation which is continued 
until a user specified tolerance is reached. In theory and in practice it offers much faster 
convergence and smaller errors than the standard trapezoidal or Simpson technique. 
[Ref. 16] 

If none of these conditions is satisfied the program defaults to a computation based 
on a numerical inverse Laplace transform. The integrand of (59) is expressed in the 


Laplace domain as [Ref. 10, equation 6.643.4}] 


Fs) =1-(54.) <exp{-(F)(1--)| ie? 


Then the integration may be determined from the inverse transform of 


{2)0-z5)}} ae 


This is useful for many intermediate values of the input arguments. 








o(a,B) 2 Fs)=4x/1-(1) xiEXp 
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D. MULTI-TONE INTERFERENCE RESULTS 


First, the results for the scenario in which at most one tone is jammed per hop slot 
are presented. The performance for a wide variety of channel conditions is computed for 
several representative signal to thermal noise ratios. All figures are for one hop per bit / 
(L=1). For each channel case there are two prime questions to be answered. First, "How 


does performance vary for fixed values of g?" Second, "Is the no fading analytic and no 





thermal noise worst case gq still valid for fading channels?" To answer the second 
question the proposed worst case value of qg, g, ,1s used to compute performance as well 
as some values of g offset from g,. To show the trend, the results of g,+ 4 and q, -4 are 
compared with g, based performance. Table 2 provides an overview of the channel 


conditions considered. 
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TABLE 2. MULTI-TONE CHANNEL CONDITIONS 


| NAME | EJN,dB |  — R,/ QT 
| Ricean Fading Ricean Fading | 










[007 Very Weak 
507 Weak 
35Ttow 
[10 7 Moderate 





507 Weak 
1 / Strong 10 / Moderate 
a J 


| Figure 35 | 20 : 10 / Moderate Tes Strong 









Initially, we consider the case where fading has very little influence to compare with 
a no fading analytic result. Figure 25 is an illustration of the performance obtained with 
a moderate signal-to-thermal noise ratio that typically provides a 10° bit error ratio in the 
absence of fading and interference. Thus, Figure 25 provides a good basis for 
comparison with the no fading case. The worst case performance is seen as the envelope 
of the fixed g curves. As the signal-to-jammer power increases, the performance quickly 
approaches 2x10”, which is very close to the no jamming performance. We can also see 
that the degradation inflicted by multi-tone jamming is strongly influenced the choice of a 


fixed number of interference tones. 
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Figure 25. E,/N,=13.35dB, R, =100, R,=100 


To bracket the range of expected performance the case of extreme of Rayleigh 
jammer fading and very weak signal fading is depicted in Figure 26. Here the 
performance is slightly improved at lower signal-to-jammer power ratios than 1s 
observed in Figure 25. Equivalently, the worst case envelope is slightly inclined more 
toward the y-axis indicating better performance. However, at low signal-to-jammer 


power levels, the performance 1s still relatively poor. 


70 


1E-2 


1E-3 


1E-4 


1E-5 





q=1 q=10 q=100 q=1000 














-30.1 -24.08 -18.06 -12.04 6.02 0.0 6.02 12.04 18.06 24.08 30.10 
Eb/NJ (dB) 


Figure 26. E,/N,=13.35dB, R,=100, R,=0.0 
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Figure 27. E,/N,=13.35dB, R,= 0, R,;=100 


At the other extreme, a pessimistic performance 1s obtained when the signal suffers 
Rayleigh fading, but the jammer has very weak fading. The jammer ts able to reduce 
performance overall and maintain a high bit error ratio (BER) for higher signal-to-jammer 
power ratios. Even at E,/N,=30dB, the effects of fading and jamming produce a BER 


near 4x10°. 
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Figure 28. E,/N,=13.35dB, R,=10, R,=10 


In Figure 28 we display the performance when both the signal and the jammer 
experience moderate Ricean fading. As expected the performance lies between the 
optimistic case of Figure 26 and the pessimistic case of Figure 27. In this case each tone 
contains the same direct-to-diffuse power ratio, but as E,/N, grows above OdB, the direct 
signal power is greater than the direct jammer power and is able to mitigate the influence 


of multi-tone jamming. 


73 


Pb 

5E-1 
q=1 q=10 q=100 q=1000 

364) eee i ie opt 








QE-1 |--\Qp--- 2-222 ee ee een ne een eee eee 


1E-1 embinc: Siig — 


SE-2 





3E-2 —- a ; 
-30.1 -24.08 -18.06 -12.04 -6.02 0.0 6.02 12.04 18.06 24.08 30.10 


Eb/NJ (dB) 


Figure 29. E,/N.=13.354B, R,=0, R,=0 


In Figure 29, the signal and jammer again share the same direct-to-diffuse power 
ratio, but now Rayleigh fading is assumed. The performance 1s worse overall compared 
with the moderate Ricean case and approximately 3dB worse at higher signal-to-jammer 
power levels. Also note that the worst case occurs at g=1000. Therefore, the greatest 
degradation occurs when the jammer is forced into the multi-tone equivalent of 


broadband jamming. 
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Figure 30. E,/N,=13.35dB, R,=10, R,=1 


In Figure 30, the channel conditions are more favorable to the communicator and 
represent a ten fold increase in the R, over the R, from the performance depicted in 
Figure 28. At the critical value of OdB for the signal-to-jammer power ratio (E,/N,) the 
performance in the favorable channel ( Figure 30) is 3.3x107 while in the moderate 
Ricean channel (Figure 28) the same E,/N, produces a BER of 2.1x10™. So the fact that 
the jammer suffers more fading over the Ricean channel produces a small improvement 


for the communicator. 


7 


q=1 q=10 g=100 q=1000 


I]E=1| Sap 





1E-2 | 


1E-3 




















1E-5 | : 
-30.1 -24.08 -18.06 -12.04 -6.02 0.0 6.02 12.04 18.06 24.08 30.10 


Eb/NJ (dB) 


Figure 31. E,/N,=13.35dB, R,=50, R,=50 


Figure 31 is similar to Figure 25, but includes less fading. The performance is 
essentially the same as the very weak fading performance, but with some slight 
improvement for E,/N, greater than OdB. As in earlier fading channel cases, the 


performance varies greatly for fixed values of g when the jammer has a power advantage. 


76 














-30.1 -15.05 0.0 15.05 30.10 
Eb/NJ (dB) 


Figure 32. E,/N,=13.35dB, R,=25, R;=25 
Figure 32 depicts another drop in the channel direct-to-diffuse ratio shared by the 


signal and jammer. At this level the difference in performance below the very weak case 
where the direct-to-diffuse ratio was 100 is noticeable. When the jammer enjoys a power 
advantage, the BER is still unacceptably high. The worst case envelope curve is nearly 


linear above E,/N,=0dB and nearly flat for all choices of g above OdB. 
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Figure 33. E,/N,=20dB, R,=10, R,=10 


In Figure 33 we begin to look at the performance for larger signal-to-thermal noise 
ratios. An increase in E,/N, from 13.35dB (Figure 28) to 20dB (Figure 33) both in 
moderate Ricean fading, provides an overall improvement in performance. However, the 
improvement is significant when E,/N, is greater than -12dB and provides nearly a factor 


of ten at OdB. 
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Figure 34. E,/N,=20dB, R,=10, R,=10 


Figure 34 shows how rapidly performance can degrade even for a strong 20dB 
signal-to-thermal noise power ratio when the signal fading progresses from moderate 
(Figure 33) to strong Ricean fading. Similar to the situation observed in Figure 27 , the 
worst case performance is obtained when the jammer places an interfering tone in every 
hop slot. The jammer therefore, need not be very sophisticated to prevent efficient 
communications. In Figure 35 the roles are reversed. With the same signal-to-thermal 
noise ratio, the signal experiences moderate Ricean fading while the jammer experiences 


strong Ricean fading. 
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Figure 35. E,/N,=20dB, R,=10, R,=1 


With moderate Ricean signal fading a particular value of g, 1s required to cause the 


worst case performance demanding greater jammer sophistication. The strongly faded 


jamming tone 1s much less effective when E,/N, is greater than OdB. Previously we 


observed an approximate tenfold drop in BER experienced for the E,/N, =13.35dB case as 


the jammer fading grew from moderate (R,=10) to strong (R,=1) at E,/N,=0dB. This 


trend is more pronounced at the E,/N, =20dB level as we transition from Figure 34 to 


Figure 35 where the improvement is greater than a factor of 25. 
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The influence of fading on the selection of g to cause worst case performance in the 
multi-tone environment remains to be answered. For FFH/BFSK in the absence of fading 
and thermal noise we expect the worst case to occur when g 1s chosen as the integer 
portion of the jammer-to-signal power ratio. A comparison of the following worst case 
performance curves will show them to be the envelope of the previous fixed g computed 


Curves. 
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Figure 36. Worst Case Multi-Tone Jamming for 
E/N, =13.35dB, R,=0. R,=100 


The most pessimistic case for the communicator is pictured in Figure 36 where the 
signal has Rayleigh fading and the jammer enjoys essentially no fading. We see that in 


this situation the jammer does not need to be very sophisticated to hinder communication. 
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The results displayed in Figure 36 and in Figure 37 indicate there is essentially no 


difference in the g, performance and the q, +4 and q, -4 cases. 
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Figure 37. Worst Case Multi-Tone Jamming 
E,/N, =13.35dB, R,=0, R,=0 


The conclusion to draw from Figure 37 is the degradation in channel conditions for 
the jammer does not significantly assist the communicator when both are acting through 


Rayleigh channels. 
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For convenience we will denote the g,+4 and gq,-4 collectively as Aq,. It is not until 
there is some appreciable direct signal power reaching the receiver that the g, and Aq, 


begin to show some deviation. 
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Figure 38. Worst Case Multi-Tone Jamming 
E,/N, =20dB, R,=10, R,=10, 


In Figure 38 the Ag, curves dip below the anticipated worst case curve. The trend 


progresses as the direct-to-diffuse ratio increases in Figure 39 and Figure 40 respectively. 
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Figure 39. Worst Case Multi-Tone Jamming 
E,/N, =13.35dB, R,=25, R,=25. 


Since the value of g is required to be a positive integer greater than or equal to one, 
there is a small discontinuity in the g,-4 curve. At those signal-to-jammer power values 


where the gp value would be less than 5, the g,-4 defaults to g=1 for the computation. 
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Figure 40. Worst Case Multi-Tone Jamming 
E,/N, =13.35dB, R,=50, R»=50, 


The salient observation in Figure 38, Figure 39, and Figure 40 is that the anticipated 
worst case performance in the no fading and no thermal noise case remains the worst case 


performance for the weakly and moderately faded channels with moderate thermal noise. 
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For comparison Figure 41, Figure 42, Figure 43 displays the composite worst case 


performance for various signal-to-thermal noise ratios. 
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Figure 41. Worst Case Multi-Tone Jamming for Several 
Values of E,/N,, R,=1, Rj=1 


Figure 41 displays the performance when both the signal and the jammer suffer 
strong fading. The diffuse nature of the channel keeps the probability of bit error 
relatively high despite high signal-to-jammer power ratios and serves more to the 


detriment to the communicator than the jammer. 
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When the channel has moderate vice strong Ricean fading the performance is that 
shown in Figure 42. Above -6dB E,/N, there is a significant improvement in 


performance over that in Figure 41. The noticeable change in slope for the 20dB curve is 
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Figure 42. Worst Case Multi-Tone Jamming for Several 
Values of E,/N,, R,=10, R,=10 


a consequence of the signal-to-jammer phase term in the bit error analysis. When the 
Signal amplitude is consistently greater than the jammer's amplitude, the phase difference 
has little effect and the rate of improving performance increases. This slope change is 


also noticeable in Figure 43 for the 20 dB curve. 
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Figure 43. Worst Case Multi-Tone Jamming for Several 
Values of E,/N,, R,=100, R»=1 


Figure 43 is a composite of performance under channel conditions favorable to the 
communicator. However, the 10:1 ratio in the signal's direct-to-diffuse over the jammer's 
direct-to-diffuse ratio provides little relief to the communicator. 

Figure 44 displays the performance when the jammer's intelligence is relaxed to 
allow two interference tones in a hop slot when neither the signal nor the jammer 


experience fading. 
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Figure 44. Worst Case Multi-Tone Interference when 
Allowing Two Jamming Tones per Hop Slot 


Figure 44 exhibits the same linear improvement in performance as E,/N, decreases 
for several values of E,/N,. There is a marked change in slope at E,/N, =0dB. Overall 
performance is better than the perfromance when at most one tone per hop slot is jammed 


in a near Gaussian channel as depicted in Figure 40. 
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IV. CONCLUSIONS 


A. FFH/MFSK PARTIAL-BAND NOISE JAMMING 


A performance analysis has been completed for the conventional fast 
frequency-hopped M-ary orthogonal frequency-shift keyed receiver employing diversity 
with noncoherent reception of signals transmitted over Ricean fading channels. We can 
draw several conclusions based on the results. 

Diversity alone is insufficient to overcome the effects of worst case partial band 
jamming acting in Ricean fading channels. The intelligent jammer can optimize the 
fraction of bandwidth jammed to force a worst case bit error ratio (BER). The 
noncoherent combining losses experienced by noncoherent systems is aggravated by the 
low signal-to-total noise ratio. These two factors combine to mitigate the advantages of 
increased diversity in a constant energy per hop system. 

Increasing the modulation order provides only modest improvements in performance 
which decrease as diversity increases. For most channel conditions, the performance for 
M=4 and M=$8 is essentially the same for all Z greater than four. The only remaining 
advantage for increased modulation order is to offset the reduction in data rate that results 
from increasing diversity, L, in a constant energy per hop system. 

When the signal experiences Rayleigh or strong Ricean fading, the jammer need not 
be very sophisticated to cause marked impairment to the communicator. Additionally, 


when comparing the partial-band jamming strategy for M=2 and Z=1 in which the signal 
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suffers moderate to strong fading to multi-tone jamming under the same conditions, we 
see that the partial-band jamming strategy causes a higher BER while requiring less 


precise knowledge of the receiver by the jammer. 
B. FFH/BFSK MULTI-TONE JAMMING 


In this thesis, the case of Z=1 is considered for the fast frequency-hopped receiver 
described above. The influence of Ricean fading of both the signal and the jammer is 
included in the analysis as is the effect of thermal noise. 

First, the counter-intuitive result is obtained that Ricean fading of the jammer 
provides little, if any, relief to the communicator, while, as expected, moderate to strong 
fading of the signal tone markedly reduces performance. Worst case multi-tone jamming 
combined with signal fading produces unacceptably high probabilities of bit error, i.e., 
greater than 10°. Above this level it is unlikely that forward error correction coding can 
be successfully implemented. 

The intelligent jammer can select the number of tones to jam which produces a worst 
case performance. The value of g chosen in the fading channel with thermal noise to 
cause worst case performance is the same value of g arrived at analytically for the noise 
free case without fading. For all combinations of channel conditions, not until the 
Signal-to-jammer power ratio exceeds 10dB is the performance essentially thermal noise 
limited. When the jammer has a 3dB or better power advantage, the performance, in even 
the most optimistic channel conditions, is very poor and is dominated by the jammer 


power. 
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Relaxing the jammer intelligence to allow the case of possibly two interference tones 
in a hop slot 1s less efficient from the jammer's point of view. Performance is worse 
when the jammer is intelligent enough to jam not more than one signal tone per hop. 
Moreover, to be effective, a multi-tone jammer needs to be more intelligent than a partial 
band jammer since the former needs to know the exact location of the signal frequencies 


in addition to the bandwidth knowledge required of the latter. 
C. RECOMMENDATIONS 


The added complexity of a weighted diversity summer (noise normalized receiver) is 
warranted in most anti-jam systems. Efforts to reduce the uncoded BER below 10° will 
likely produce better overall performance when forward error correcting codes are 
implemented. When available, side information should be used to avoid a jammer or 
discard jammed hops. 

The partial-band noise jammer who 1s intelligent enough to force the worst case 
scenario is more threatening to the communication system than multi-tone interference 
and will require a detailed threat analysis for hostile sources. Though "dumb", 
non-hostile interference sources also warrant attention in mission planning since even a 


single signal tone consistently jammed can cause detriment to the communicator. 
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